Comparison of CORSIKA and COSMOS simulations 
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Abstract. Ultra-high-energy cosmic rays (UHECRs) refer to cosmic rays with energy above 10 18 eV. UHECR experiments 
utilize simulations of extensive air shower to estimate the properties of UHECRs. The Telescope Array (TA) experiment 
employs the Monte Carlo codes of CORSIKA and COSMOS to obtain EAS simulations. In this paper, we compare the results 
of the simulations obtained from CORSIKA and COSMOS and report differences between them in terms of the longitudinal 
distribution, Xmax-value, calorimetric energy, and energy spectrum at ground. 
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I. INTRODUCTION 



II. SIMULATIONS 



Ultra-high-energy cosmic rays (UHECRs) with energy 
greater than 10 18 eV arrive at Earth from space. Exper- 
iments to detect UHECRs utilize extensive air shower 
(EAS) to estimate their energy, composition and arrival 
direction. In an EAS, the cascade of interactions, induced 
by a UHECR in the upper atmosphere, results in a very 
large number of secondary particles: of order 10 12 parti- 
cles for a primary particle with energy 10 19 eV. 

Monte Carlo codes have been used to simulate EAS. 
The codes reproduce cascades of particles initiated by the 
interaction between primary UHECRs and atmospheric 
nuclei. EAS simulations give the spatial, temporal, en- 
ergy, and angular distributions of secondary, air shower 
particles. To extract the information of primary particles 
in experiments, that is, the energy, composition and ar- 
rival direction, the measured quantities of EAS need to be 
compared with those from simulations. Hence, the accu- 
rate reproduction of EAS is an essential part of UHECR 
experiments. 

CORSIKA [1] and COSMOS [2] are among the Monte 
Carlo codes. Most experiments including the recent one 
at the Pierre Auger Cosmic Ray Observatory have re- 
lied simulations from CORSIKA. On the other hand, the 
Telescope Array (TA) experiment [3] has employed both 
CORSIKA and COSMOS, to cross-check the simula- 
tions from the codes. In this paper, we compare COR- 
SIKA and COSMOS simulations by examining the dif- 
ferences in the quantities such as the longitudinal dis- 
tribution, Xmax-value, calorimetric energy, and energy 
spectrum at ground in EAS. 



We employed the most recent versions of codes: Ver- 
sion 7.54 for COSMOS and Version 6960 for CORSIKA. 
Each code has an option to set interaction models; dif- 
ferent interaction models result in somewhat different 
results. We chose the following interaction models: 1) 
for CORSIKA, QGSJETII-03 [4] for high-energy (above 
80 GeV) hadronic interactions, FLUKA (v.2008.3c) for 
low-energy (below 80 GeV) hadronic interactions, and 
EGS4 for electromagnetic interactions, 2) for COSMOS, 
QGSJETII-03 for high-energy (above 80 GeV) hadronic 
interactions, PHITS [5] and JAM [6] for low-energy (be- 
low 80 GeV) hadronic interactions, and Tasi's and Nel- 
son's formula for electromagnetic interactions. Please 
see references for CORSIKA and COSMOS for details. 

The above interaction models are the default, so most- 
widely used models. In this paper, we intend to compare 
the results of most-commonly employed CORSIKA and 
COSMOS simulations. We note that even with the same 
interaction models, the results of CORSIKA and COS- 
MOS simulations could be different, because of the dif- 
ferences in handling the development of EAS. Compar- 
isons of CORSIKA and COSMOS simulations with dif- 
ferent interaction models will be reported in an upcoming 
journal paper. 

We adopted exactly same parameters for CORSIKA 
and COSMOS simulations. For thinning parameter, 10~ 7 
was used (Hillas thinning algorithm) [7]. The ground 
level was located at 875 g/cm 2 (1430 m), and the atmo- 
spheric depth and the Earth magnetic field suitable for 
the TA site were applied. 
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FIGURE 1. Longitudinal development of proton EAS with 
the primary energy 10 I9 eV and the zenith angle degree. 
Panels show the numbers of photons, electrons, muons, and 
hadrons as a function of atmospheric depth. Black lines are the 
CORSIKA results, while blue lines are the COSMOS results. 
The lines are averages of 50 shower events. 



We generated air shower events with cosine of the 
zenith angle of 1, 0.95, 0.9, • • •, 0.7 and primary energy 
of 10 18 ' 5 , 10 18 - 75 , 10 19 , 10 20 ' 25 eV for proton and iron 
primaries. Altogether, about 10,000 showers were gener- 
ated with each of CORSIKA and COSMOS. The results 
below are based on some of the shower simulations. 



III. COMPARISONS 

1. LONGITUDINAL DISTRIBUTION 

Figure 1 shows a typical longitudinal development of 
air showers obtained with CORSIKA and COSMOS. 
Overall, the numbers of secondary particles are predicted 
to be larger with CORSIKA than with COSMOS. There 
are noticeable differences. For instance, the maximum 
difference in the photon number reaches up to ~ 10%. 
The difference in the hadron number is much larger, 
although the number itself is much smaller than those of 
other particles. On the other hand, the agreement in the 
electron and muon numbers is in a acceptable level; the 
maximum difference in the electron number is ~ 5% or 
less. 



Vertical shower 
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FIGURE 2. Xmax as a function of the primary energy. The 
results shown are for vertical showers with zenith angle 
degree. The upper part is for proton primaries, and the lower 
part is for iron primaries. Black dots are the CORSIKA results, 
while blue dots are the COSMOS results. Those are averages 
of 50 shower events. Lines are least chi-square fits of points. 



2. Xmax 

Xmax is the atmospheric depth of the shower maxi- 
mum, specifically the maximum of the electron distribu- 
tion. It is one of most important quantities in UHECR 
experiments; it is mainly used to determine the com- 
position. We employed the Geisser-Hillas function [8] 
to determine the longitudinal development and Xmax, 
which has been widely employed in other studies. Fig- 
ure 2 shows Xmax for both proton and iron primaries 
from COSMOS and the CORISKA simulations. In gen- 
eral, iron-initiated air showers penetrate less than proton- 
initiated air showers, so Xmax for iron primaries is less 
than that for proton primaries. And Xmax increases with 
the energy of primaries. The agreement in COSMOS and 
the CORSIKA simulations is reasonably good; the differ- 
ence in Xmax is only a few percent at most. 



3. CALORIMETRIC ENERGY 

As EASs develops, a part of the energy of primary 
particles (Eq) is deposited into air molecules and even- 
tually radiated as fluorescence lights. But a fraction of 
the energy is carried away by secondary particles, not 
contributing to fluorescence lights. A correction for the 
so-called missing energy must be applied to the mea- 
surement of the calorimetric energy (£«,/), in order to 
correctly determine the primary energy, Eq, from obser- 
vation of fluorescence lights. Here, we calculated the 
missing energy by following the prescription described 
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FIGURE 3. Calorimetric energy as a function of the primary 
energy for proton primaries. The results shown are for vertical 
showers with zenith angle degree. Black dots are the COR- 
SIKA results, while blue dots are the COSMOS results. Those 
are averages of 50 shower events. Lines are least chi-square fits 
of points. 



by [9], and so the calorimetric energy. Figure 3 shows 
the resulting E ca \ for proton primaries from CORSIKA 
and COSMOS simulations. Our result indicates that E ca i 
from COSMOS is ~ 2% larger than that from COR- 
SIKA. This implies that the primary energy estimated 
with fluorescence detector would be ~ 2% smaller when 
COSMOS simulations are used. 
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FIGURE 4. Energy (rest-mass energy no included) distribu- 
tion of photons at ground over the entire ground for proton EAS 
with the primary energy 10 19 eV and the zenith angle degree. 
Black lines are the CORSIKA results, while blue lines are the 
COSMOS results. The lines are averages of 50 shower events. 




4. ENERGY DISTRIBUTION AT THE GROUND 

As a consequence of EAS, a number of secondary par- 
ticles arrive at ground. In experiments, those particles are 
registered by surface detectors and used to estimate the 
energy and arrival direction of the primary particles. In 
Figures 4, 5, and 6, we compare the energy (rest-mass en- 
ergy no included) distribution of photons, electrons, and 
muons reached at ground from CORSIKA and COSMOS 
simulations. Particles in the core are included. The figure 
indicates a good agreement between the CORSIKA and 
COSMOS results; the typical difference is ~ 3% or so. 
Although not shown here, the difference in the distribu- 
tion of hadrons is much larger. But again, the number and 
energy of hadrons are much smaller than those of other 
particles. 

IV. DISCUSSION 

Monte Carlo codes of CORSIKA and COSMOS are cur- 
rently used to analyze the data of the TA experiment. In 
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FIGURE 5. Energy (rest-mass energy no included) distribu- 
tion of electrons at ground over the entire ground for proton 
EAS with the primary energy 10 19 eV and the zenith angle 
degree. Black lines are the CORSIKA results, while blue lines 
are the COSMOS results. The lines are averages of 50 shower 
events. 



this paper, we compared simulations of EAS using COR- 
SIKA and COSMOS codes and quantified the differences 
in the simulations. For the longitudinal distribution of 
photons, electrons, and muons, we found the maximum 
difference of ~ 10%. For Xmax and calorimetric en- 
ergy, the difference is typically a few percent. The differ- 
ence in the energy distribution of photons, electrons, and 
muons at ground is again typically a few percent. This 
implies that we should expect an uncertainty of a few per- 
cent in the estimation of the primary energy in UHECR 
experiments including the TA experiment, solely due to 
the uncertainty in EAS simulations. 

Finally, it is worthwhile to mention that there is a 
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FIGURE 6. Energy (rest-mass energy no included) distribu- 
tion of muons at ground over the entire ground for proton EAS 
with the primary energy 10 19 eV and the zenith angle degree. 
Black lines are the CORSIKA results, while blue lines are the 
COSMOS results. The lines are averages of 50 shower events. 



large difference in the production of secondary hadrons 
in the CORSIKA and COSMOS codes. We found that the 
production is rather sensitive to the low-energy (below 80 
GeV) hadronic interaction model. This may indicate that 
low-energy hadronic interaction models in Monte Carlo 
codes for EAS need to be further investigated. 
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